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Abstract 

A schematic model is presented that allows one to study the behavior of interacting pseudo- 
Fermi matter, locked in a thermostatic box. As a function of the box volume and temperature, 
the matter is seen to show all of the familiar characteristics of a Van der Waals gas, which 
include the coexistence of two phases under certain circumstances and the presence of a critical 
point. 

PACS numbers: 21.65+f, 21.60.Ev, 25.70.Pq 
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I. INTRODUCTION 

The possibility that a hquid-gas phase transition in finite nuclei may manifest it- 
self via copious production of intermediate-mass fragments (IMF) in energetic heavy-ion 
reactions has driven both, theoretical and experimental studies of nuclear multifragmen- 
tation over more than a decade. A prominent role in theoretical considerations is played 
by the concept of a freezeout configuration. |jl|, ^ This concept implies the existence of a 
definite volume, within which the system reaches a state close to thermal equilibrium. 
While the existence of an effective freezout volume may be debatable, the concept itself 
is useful for understanding the possible behavior of nuclear matter under various con- 
ditions. With such a didactic strategy in mind, and as an extension of earlier studies 
modelling the behavior of finite Fermi systems, |P, |, |^ the present study considers 
pseudo-Fermi matter confined to a box of definite volume. It evaluates the isothermal 
behavior of such matter and its dependence on box volume and temperature. The ut- 
most simplicity of the formalism allows one to gain insight into physical phenomena that 
may be obscured in more rigorous approaches. 

II. THEORETICAL FORMALISM 

The present study considers a scenario of nuclear pseudo-Fermi matter of mass num- 
ber A, locked in a spherical box of a volume V and kept at constant temperature T. 
In the proposed formalism, thermostatic properties of the matter are modelled by two 
equations, the isochoric caloric equation of state, and the zero-temperature equation of 
state. The equilibrium state of the system is then found as a function of V and T, based 
on the requirement that the free energy of the system be minimal. 

First, we consider homogeneous systems with constant density throughout their vol- 
ume, i.e., systems in states where only one phase is present. The isochoric caloric equa- 
tion of state for such a system is taken in a simple form adequate for low-temperature 
Fermi gases 
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where E^f^^^^ is the thermal energy and a is the level density parameter. The latter 
parameter is assumed to depend on matter density as 

a = ao{-)-h (2) 
po 

in accordance with the low-temperature Fermi gas model. 

It is worthwhile keeping in mind, that for Fermi gases, Eq. |1| is a good approximation 
only for temperatures that are small compared to the Fermi energy, i.e., for T << Ep- 
Notably, for diluted or very hot Fermi systems, with the matter density approaching zero, 
p— > 0, or for high temperatures, T— > oo, the caloric equation of state approaches 
asymptotically that of a classical gas: 

3 

E therm = 2^^"^ 

For the sake of simplicity and without loss of generality, the present study uses Eq. 
over the full range of matter densities and temperatures considered. The term "pseudo- 
Fermi" matter is used to distinguish the matter considered here from true Fermi matter. 

The second defining equation, the zero-temperature equation of state, expresses the 
compressional (potential) energy of the system as a function of the system volume or as 
a function of matter density. The present study adopts the harmonic approximation in 
the form used in the Expanding Emitting Source Model (EESM) P], whereby the in- 
medium nucleonic (potential) energy changes quadratically with the relative deviation 
of the actual matter density from the ground-state density, by an amount 

^compr = — £^(1 —)'^. (4) 

Po 

In Eq. ecompr and €3 are compressional and ground-state binding energies per nu- 
cleon, respectively, and p and po are the actual and the ground-state matter densities, 
respectively. Equation ^ implies an effective ground-state incompressibility constant of 
Ko=-l8eB- Assuming eB=-8 MeV, the effective incompressibility, including the effects 
of surface tension, equals Ko = 144 MeV. Note, that within the harmonic approxima- 
tion, for infinite nuclear matter characterized by ~ —16 MeV, the incompressibility 
constant is Kg ~ 288 MeV. The latter value places the present harmonic approximation 
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"neutrally" between the currently considered limits of "soft" and "hard" equations of 
state for nuclear matter. 

While the two defining equations |1| and ^ may be considered rather crude approxima- 
tions, they do contain the essential physics responsible for first-order phase transitions 
and critical phenomena. Given these two equation, one can write expressions for all 
thermodynamic quantities characterizing the system, including the Helmholtz free en- 
ergy F. The state of the system can then be found by minimizing the free energy, for 
any box volume V and temperature T. 

Based on Eq. |l], one can write for the entropy 5 of a homogeneous system 

S= l.de = 2\^. (5) 
Jo T 

The free energy, F, for a homogeneous, single-phase system is given by 

^ ~ ^total '~ ~ ^compr + ^therm ~ ^^T" = -Ecompr ~ • (6) 

In Eq. ^, for sake of simplicity, the free energy is given relative to a constant ground 
state energy of the system. Further, using Eqs. |^ and ^, Eq. ^ can be rewritten in a 
form revealing explicitly the important dependence on the matter density p, i.e., 

F = En{l--^r-a^{^r'^T'. (7) 

Po Po 

Eq. 1^ allows one to write equations for the system pressure p and the chemical po- 
tential /i as functions of volume (matter density) and temperature. 

The pressure p can be expressed generally as the negative partial derivative of the free 
energy F with respect to volume V, at fixed number of nucleons A and fixed temperature 
T, i.e., 



Thus, for the case of pseudo-Fermi matter considered here, one obtains based on 
Eqs. and | 

p = 2esPo{l - ^){^f + \aoPo{^f^T\ (9) 

Po Po 3 Po 

where and ao are the binding energy per nucleon and the level-density parameter per 
nucleon, respectively. 
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The chemical potential /i can be expressed generally as a partial derivative of free 
energy F with respect to the number of nucleons A, taken at fixed volume V and fixed 
temperature T, i.e., 

OF 

Using Eqs. ^ and |10|, and noting further that p = A/V , Eb = Aes and Oo = Aao-, 
one obtains for the chemical potential 

^^ = eB[l-^^ + K^f]-\a^{^)-iT\ (11) 

Po Po 3 Po 

It is worth noting in Eq. |Tl| that, as a result of the requirement for V to be constant, 
the magnitude of the chemical potential differs significantly from the value of the free 
energy per nucleon. For example, the contribution of thermal excitation to the chemical 
potential is only one-third of what constitutes the thermal part of the free energy per 
nucleon. 



III. LIQUID-GAS COEXISTENCE 

When confined to a thermostatic box, matter will eventually fill all of the available 
volume V so as to minimize the free energy of the system. For non-interacting matter, 
the minimum free energy always corresponds to a uniform density distribution. This is 
generally not true for an interacting system. In particular, for interacting pseudo-Fermi 
matter at the minimum free energy, low-density and high-density phases coexist. 

For a two-phase system of A nucleons confined to volume V at temperature T, the 
free energy can be expressed as a function of two variables, e.g., in terms of the volume 
of the gaseous phase, Vgas, and the number of nucleons contained in this phase, Agas- 

F -Fgas ~l~ Fliquidi (■^^) 

Inserting for Fgas and Fuguid the expressions given by Eq. ^ for the corresponding 
numbers of nucleons Agas and Auguid, respectively, one obtains 

Fgas = AgaM^ - - Agasa^i^r'^T' (13) 

'gasPo ^gasPo 
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FIG. 1: Isotherms calculated for an ^^''Au-like system. Dotted lines illustrate isotherms for 
hypothetical single-phase matter, while the dashed line visualizes the boundary of the liquid- 
gas coexistence domain. 

and 

The condition for the minimum free energy can be expressed in form of two equations, 
reflecting requirements of dynamical and chemical equilibrium, respectively 



)T,Aga, — )T,Aaa, ~ y^y. ~)T,Auquid — Pgas — Pliquid — U l^iOJ 



gas ^'gas ^''liquid 



and 



' 9F OFgas fdFiiquids „ 

.^1 )T,Vgas - [tTI )T,Vgas - )T,Vug^,a - l^gas - ^^liquid - U. {IK)} 

(^-^gas ^-^gas ^^-^liquid 

Results of a numerical minimization of the free energy of a two-phase system are 
shown in Figs. Fig. 1 presents system isotherms, as predicted by the present for- 
malism, in a representation of system pressure p versus system volume V. Note that 



Eqs. |T3| and II represent single-phase states as special cases. A pure liquid/gas state is 




FIG. 2: Liquid-gas coexistence line in the temperature versus matter density representation, 
as predicted by the present formahsm. 

thus among the possible outcomes of the calculations with Agas/uquid = 0. As seen in 
Fig. D, the isotherms feature segments representing both, pure liquid and pure gas and, 
notably, the liquid-gas coexistence "plateaus". The presence of coexistence plateaus 
does not come as a surprise as the harmonic interaction term of Eq. ^ has the salient 
characteristics of a Van der Waals interaction. 

It is worth noting that the coexistence plateaus in this calculation result naturally 
from the actual minimization of the free energy and are not obtained via the well- 
known phenomenological Maxwell construct. Isotherms for hypothetical single-phase 
states are shown as dotted lines. These latter isotherms feature domains of spinodal 
instability characterized by negative compressibility. The dashed line in Fig. |l| illustrates 
the boundary of the liquid-gas coexistence domain. The "summit" point of this curve 
represents the critical point for the system and corresponds to a critical temperature 
of Tc ~ 10.0 MeV. At temperatures higher than Tc, the system can reside only in a 
single-phase, vapor state. 

A different representation of the liquid-gas coexistence line is illustrated in Fig. |^. In 
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FIG. 3: Isochoric caloric curves for the boxed pseudo-Fermi matter, calculated for different 
volumes of the confining box, as indicated by labels. 

this case, the temperature T is plotted versus the matter density, for points along the 
boundary of the liquid-gas coexistence line (dashed line in Fig. |I|). In the domain below 
this curve, the system is in a two-phase state. In this coexistence domain, the densities 
of gaseous and liquid phases at a given temperature are depicted by the intersection 
points of the coexistence curve with the horizontal line T=const. On the left shoulder 
of the curve and in the domain further left to it, the system is in a pure gaseous state, 
while on the right shoulder and in the domain further right to it, the system is in a pure 
liquid state. The difference between hquid and gaseous phases vanishes for temperatures 
above T = T^ = 10.0 MeV. 

The search for the minimum free energy for a given temperature T and given volume 
V yields, along with the volumes Vgas/iiquid and densities Pgas/iiquid of gaseous and liquid 
phases, the total excitation energy E*^^^i of the boxed matter. This allows one to con- 
struct the isochoric caloric curves for the modelled system at different volumes of the 
confining box. A set of such caloric curves is illustrated in Fig. ^. As expected, these 
curves do not feature plateaus reported in various experiments, but rather inconspic- 
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uous kinks at the locations on the boundary of the coexistence domain. 



IV. CRITICAL BEHAVIOR 



One of the sahent features of critical behavior in Van der Waals systems is the presence 
of a singularity at the critical point. Here, when the system temperature T approaches 
the critical temperature from below, the difference between the densities of coexisting 
liquid and gaseous phases vanishes according to a power law 

Pliquid — Pgas = C{Tc — T)^ , (17) 



In Eq. C is a constant and (3 is the critical exponent. The magnitude of the critical 
exponent can be extracted conveniently by fitting a straight line to a double-logarithmic 
plot of {piiquid — Pgas) vcTsus (T^ — T). Such a plot obtained in present calculations is 
shown in Fig. ^ The plot features, indeed, an approximately 2-MeV wide linear domain 
extending from the critical temperature down to lower temperatures. A linear fit 
to this domain allows one to extract the "coordinates" of the critical point of Tc= 10.0 
MeV, and pc/po=0. 42. Furthermore, it yields a value of /3=0.51 for the critical exponent. 

Given the very schematic nature of the present formalism, the predicted characteris- 
tics of the critical point are well within a reasonable domain that can be inferred from 
more sophisticated, but less transparent calculations. 



V. DISCUSSION 



A simple formalism has been presented that allows one to model the behavior of inter- 
acting pseudo-Fermi matter under the conditions of controlled volume and temperature. 
The formalism is shown to capture the essential physics of a first-order liquid-gas phase 
transition. The calculation demonstrates the characteristics of nuclear liquid-gas coex- 
istence in a certain domain of system parameters and the presence of a critical point. 
While such characteristics are well expected, based on the the similarity of the utilized 
form of interaction to that of the Van der Waals interaction, the present formalism of- 
fers many didactical benefits. For example, due to its simplicity it offers a clear insight 
into physical phenomena it purports to describe and a relatively strict test-bench for 
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FIG. 4: Critical behavior of densities of coexisting liquid and gaseous phases at temperatures 
close to critical temperature Tg. 

possible qualitative or "hand-waving" arguments. The results obtained, while almost 
trivial, may alert one to possible challenges found by more complete models. As an ex- 
ample, it is evident from Fig. 1 that models of nuclear multifragmentation such as SMM 
and MMMC tend to overpredict fragment multiplicities significantly. This is so, 
because these models assume that the liquid phase is always at ground-state density 
Po- This forces the pressure of the surrounding gas within the freezeout volume to be 
higher than would be necessary to counter the lower pressure of the expanded liquid 
at equilibrium density. A higher gas pressure results in an overestimate of particle and 
fragment multiplicities. 

The formalism presented here leaves ample room for further refinements such as a 
more strict modelling of the diluted Fermi matter, or the inclusion of isotopic effects. At 
any rate, it offers a convenient didactic tool to achieve a better understanding of nuclear 
thermodynamics. 
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